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ABSTRACT 

We have assembled a homogeneous database of 417 ultraviolet (UV) ex- 
tinction curves for reddened sightlines having International Ultraviolet Explorer 
(lUE) spectra. We have combined these with optical and 2MASS photometry 
allowing estimates of the ratio of total-to-selective extinction, Ry, for the entire 
sample. Fitzpatrick-Massa (FM) parameters have also been found for the entire 
sample. This is the largest study of parameterized UV extinction curves yet pub- 
lished and it covers a wide range of environments, from dense molecular clouds 
to the diffuse interstellar medium (ISM), with extinctions Ay ranging from 0.50 
to 4.80. It is the first to extend far beyond the solar neighborhood and into the 
Galaxy at large, with 30 sightlines having distances > 5 kpc. Previously, the 
longest sightlines with FM parameters and Ry extended ~ 1 kpc. We find that 
(1.) the CCM extinction law applies for 93% of the sightlines, implying that dust 
processing in the Galaxy is efficient and systematic; (2.) the central wavelength 
of the 2175 A bump is constant; (3.) the 2175 A bump width is dependent on 
environment. Only four sightlines show systematic deviations from CCM, HD 
29647, 62542, 204827, and 210121. These sightlines all sample dense, molecule- 
rich clouds. The new extinction curves and values of Ry allow us to revise the 
CCM law. 

Subject headings: dust, extinction 



1. Introduction 

Astronomical ultraviolet studies were not possible before the advent of rocket-borne 
instruments and sateUites due to atmospheric absorption and Rayleigh scattering. The first 
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satellite capable of UV observations, OAO-2, was launched in 1968, and was followed by 
a host of other satellites (Cox 2000), including Thor Delta-1 (TD-1), Astronomy Nether- 
lands Satellite (ANS), and lUE, among others. The data obtained by these instruments 
are responsible for the vast majority of contributions to studies of interstellar dust, as dust 
grains preferentially extinguish short-wavelength light. By comparing heavily extinguished 
stars against their intrinsic fluxes (the pair method), we can better understand grain prop- 
erties and composition. This is important for two reasons. First, solid knowledge of grain 
properties will allow us to build more realistic grain models and correct for the effects of 
dust in stellar and galactic photometry and spectra more accurately. Second, dust grains 
are extremely sensitive to their environment, and can yield much information about local 
conditions and the physical and chemical processes which affect grains. 

The pair method, first used by Bless & Savage 1970, has since been used by many others 
to make extinction curves and study the ISM in the Milky Way (e.g., York et al. 1973; Nandy 
et al. 1975, 1976; Koorneef 1978; Seaton 1979; Snow & Scab 1980; Kester 1981; Meyer & 
Savage 1981; Aiello et al. 1982; Massa, Savage & Fitzpatrick 1983; Massa & Savage 1984; 
Savage et al. 1985; Fitzpatrick & Massa 1986, 1988, 1990, hereafter FM86, FM88, FM90; 
Clayton & Fitzpatrick 1987; Aicllo ct al. 1988; Cardelli et al. 1988, 1989, hereafter CCM; 
Papaj, Krelowski, & Wegner 1991; Papaj & Krelowski 1992, and Jenniskens & Greenberg 
1993, hereafter JG93) as well as in other galaxies (Borgman et al. 1975; Nandy et al. 1981; 
Koorneef & Code 1981; Clayton & Martin 1985; Fitzpatrick 1985, 1986; Clayton et al. 1996; 
Bianchi et al. 1996; Gordon & Clayton 1998; Misselt et al. 1999; Gordon et al. 2003). 

There are significant differences between dust grain properties in the Galaxy and those 
found in the Small and Large Magellanic Clouds (SMC and LMC, respectively). The SMC 
Bar has very strong, linear FUV extinction, while the sighthnes toward the LMC2 Supershell 
and the average LMC curve have a weaker bump and stronger FUV extinction than the 
average Galactic Ry=3.1 curve. There are also significant differences in metallicities and the 
gas-to-dust ratios between the Galaxy, LMC, and SMC (Bohlin et al. 1978; Gordon et al. 
2003; Luck & Lambert 1992). 

A major leap forward in understanding interstellar extinction occurred in 1989, when 

CCM suggested a mean extinction relation that depended on one parameter, the ratio of 
total-to-selective extinction, Ry (=Ay/E(B-V)). With only very few exceptions. Galactic 
extinction curves tend to follow this law within the uncertainties of the calculated Ky values 
and the extinction curves (Fitzpatrick 1999; Clayton et al. 2000). The Galactic diffuse ISM 
is well described by a curve where Ry = 3.1. Ry can also be thought of as a rough indicator 
of grain size, with low Ry sightlines having more small grains than high Ry sightlines. The 
CCM extinction relation generally does not apply outside the Galaxy, although recent work 
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has shown that there are sighthnes toward the LMC which may follow CCM (Gordon et 
al. 2003). This law essentially replaced the Galactic average extinction curve (Seaton 1979) 
with a family of Ry-dependent extinction curves, with each curve representing a mean curve 
for sightlines of that particular value of Ry. This also showed that many of the so-called 
"anomalous" sightlines (e.g. Clayton & Fitzpatrick 1987) were in fact normal. In addition 
to finding this relationship, CCM pointed out the usefulness of normalizing extinction curves 
to Ay, rather than E(B-V), as was usually done. Compared to E(B-V), Ay is a more basic 
quantity, and it is a direct measure of optical depth along a line of sight. Galactic extinction 
curves can be fit with six parameters (Fitzpatrick & Massa 1990) which allow for quantitative 
analysis of extinction characteristics. 

Despite all the work with extinction curves that has been done, a homogeneous database 
of Milky Way extinction curves does not exist. The many previous studies, hsted above, used 
data from different instruments, and variations of the pair method to construct extinction 
curves. So, the results of these studies cannot be compared easily. Also, these sightlines 
make up only about half of reddened sightlines available in the lUE Archive. To illustrate 
the disparateness of previous works. Table 1 compares various aspects of previous studies to 
the present work. Column (1) lists authorship, (2) lists the number of sightlines, (3) lists the 
instrument the data was from, (4) shows if FM parameters were found, (5) shows whether 
the extinction curves were constructed using the MK classification or UV stellar features, 
and (6) shows if IR photometry, specifically Ry, was considered for the entire sample. 

In order to solve this problem, we have built a homogeneous database of over four 
hundred extinction curves constructed in a uniform manner. Over 150 of these sightlines 
extend beyond 2 kpc, and 30 extend beyond 5 kpc, thus sampling a much larger volume 
of the Galaxy than ever before. This is a much larger and more complete database than 
any done previously. Until the recent release of the 2MASS database, JHK photometry was 
available for only a small fraction of the stars in the sample. Now, estimates of Ry are 
available for the entire sample. Thus, it is now possible to draw a more coherent picture of 
extinction in the Galaxy. In the following sections, we describe our database of extinction 
properties, present the results of our analyses, and discuss the implications for dust grains. 

2. UV Data and Extinction Curves 

For the sake of homogeneity, Hiltner's UBV photometry (Hiltner 1956; Hiltner & John- 
son 1956) was preferred, when possible, and was obtained for 192 sightlines. In his massive 
compilation and assessment, Nicolet (1978) found that Hiltner's photometry was consistently 
of high quality. Systematic errors were typically less than 0.01 mag (Hiltner 1956) and pho- 
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tometry agreed well with the works of others (Nicolet 1978). Photometry in the JHK bands 
was available for all the stars in the sample from the 2MASS database (Cutri et al. 2003). 

Reddened stars, observed by lUE, having spectral types 03 to B5 were selected for 
the database. An exception is HD 29647, a B8 III star, included because it has been well 
studied in the past (e.g., CardeUi & Savage 1988). This range of spectral types was chosen 
to minimize the effects of spectral mismatch in their extinction curves (Massa et al. 1983). 
Thorough discussions of the uncertainties in the pair method can be found in Massa et al. 
(1983), Savage et al. (1985), Aiello et al. (1988), and Gordon & Clayton (1998). Only 
"normal" stars were included which had good unreddened UV comparison stars. The lower 
limit on E(B-V) was 0.20. 

lUE spectra for all of the stars in the sample were obtained from the Multimission 
Archive at Space Telescope (MAST). The archive spectra were reduced using NEWSIPS 
and then recalibrated using the method developed by Massa & Fitzpatrick (2000). The 
signal-to-noise of the NEWSIPS lUE spectra have been improved by 10-50% over those of 
lUESIPS lUE spectra (Nichols & Linsky 1996). Low dispersion LWR/LWP and SWP spectra 
were selected, from either aperture. Multiple spectra from one camera were averaged and 
then the long and short- wavelength segments were merged at the shortest wavelength of the 
SWP. The wavelength coverage is ~1200 - 3200 A with a resolution of ~5 A. 

The standard pair method, in which a reddened star is compared with an unreddened 
one of the same spectral type, was used to construct each sightline's extinction curve (Massa 
et al. 1983). The comparison stars were selected from CardeUi, Sembach, & Mathis (1992) 

and dereddencd. The spectral matches were made on the basis of comparing the UV spectra 
of pairs of stars rather than matching their visible spectral types. Ry was estimated from 
the JHK colors as described in Fitzpatrick (1999). Ky was found using Ry and E(B-V). In 
general, our UV spectral classifications are in good agreement with those which have been 
published previously (e.g., Aiello et al. 1988; Clayton & Fitzpatrick 1987; Papaj et al. 1991). 
Photometry, sources, MK and UV spectral types, reddenings, and calculated Ry values for 
all sighthnes can be found in Tables 2, 3, and 4. 

The resulting extinction curves were fit with the Fitzpatrick-Massa (FM) parameteri- 
zation (Fitzpatrick & Massa 1990) and then normahzed to Ay. Setting x = 1/A, the FM 
fitting function is given by 



K^) = ^( ^ _ = Cl + C2X + C3D{X, 7, Xo) + C4F{x) (1) 



where 
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and for x > 5.9 /im 



-1 



F{x) = 0.5392(x - 5.9)^ + 0.05644(x - 5.9)^ (3) 
while F{x) = for x < 5.9 //m~^ 

Eqn. 1 results in 6 parameters, Ci, C2, C3, C4, xq, and 7, each of which describe different 
attributes of the curves. The first two, Ci and C2, account for the intercept and slope of 
the linear background. They are not independent of each other and can be merged into one 
parameter (Carnochan 1986, FM88), though that is not done here. The quantities C3 and 
C4 correspond to the strength of the bump and the curvature of the FUV rise. The last 
two parameters, Xq and 7, correspond to the central wavenumber and width of the bump, 
respectively; xq does not vary greatly from sightline to sightline, so it may be possible to 
reduce the number of parameters to four overall (Fitzpatrick 1999). As can be seen in Eqns. 
2 and 3, the bump is fit by a Drude profile and the FUV extinction is fit by a nonlinear 
function. The parameters were found using the three-step method of Gordon et al. (2003), 
over the wavelength range 2700 - 1250 A. The FM fit is not rehable longward of 2700 A (E. 
Fitzpatrick 2002, private communication). The spectra are cut at the blue end at 1250 A in 
order to exclude the Lya feature at 1215 A. The normalization of the FM parameters was 
converted from E(B-V) to Ay (JG93): 

A{V) Rv 
so that the FM parameters become 



cf^ = ci/Rv + 1.0 
.A 



cf^c,/Rv, J = 2,3,4 



The resulting normahzed FM parameters are listed in Table 5. The dust environment 
was assessed by considering the traditional density measure, Ay/d; JG93 found that sight- 
hnes passing through dense regions had Ay/d > 0.9 mag kpc~^, while more diffuse sightlines 
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had Ay/d < 0.9 mag kpc~^. We have adopted this criterion for this study. This parame- 
ter can be inaccurate since the density is averaged over the entire sighthne (JG93). Other 
groups have noted that trends in extinction characteristics, especiaUy those associated with 
the bump, can be hnked to environmental influences (FM86; CardeUi & Clayton 1991; JG93). 

Figure 1 shows the distribution of previously published sightlines in the Galaxy with 
FM parameters normalized to E(B-V) and the sightlines in this new database. Because of 
our much larger sample, we have more complete coverage of nearby associations and more 
distant sightlines than previous works, pushing out further into the Galaxy, and sampling a 
wider variety of environments. As distances were estimated using spectroscopic parallaxes, 
the uncertainties for stars plotted in Fig. 1 are large, with an average uncertainty of 50%. 

3. Discussion 

3.1. Continuum Extinction 

The relationship between the various FM parameters were investigated as well as their 
dependence on Ry~^. These are shown in Figs. 2-8. For pairs of parameters with significant 
correlations, best-fit lines are plotted, using the least absolute deviation method as this is 
less affected by outliers than minimization. 

There is a clear correlation between Ci/Ry + 1.0 and C2/Ry (see Table 6 and Fig. 2). 
This has been noted before (Carnochan 1986; FM86; JG93). According to FM88, Ci = -3.00 
C2 + 2.04, while JG93 find that Ci = (-3.11±0.11) C2 + (2.14±0.07). This relation shows 
how tightly constrained the linear component of the extinction is with respect to x (FM88). 
There is also a weak correlation between C2/Ry and C3/Ry. Neither JG93 nor FM88 found 
any apparent relationships between the parameters C2 and C3. Similarly, no correlation was 
found in our sample between these parameters. A comparison of the parameters found in 
this study and those of JG93 and FM88 are shown in Fig. 3. The values of C2 and C3 from 
the other groups have been divided by the values of Ry found in this work, and the best-fit 
line in the lower panel was found by considering points from both studies. It was similar 
to that found here: Cs/Ry = 1.77 C2/Ry + 0.67, versus Cs/Ry = 1.75 C2/Ry + 0.54. As 
both C2/Ry and Ca/Ry increase with Ry~^, and thus are infiuenced by similar environmental 
effects, perhaps it is not surprising that they themselves are correlated with each other. 

The parameters C4/Ry and C2/Ry do not show a correlation; neither do the parameters 
C2/Ry and 7. These agree with the findings of FM88 and JG93. There is also no correlation 
(r=0.13) between C4/Ry and 7. This is in disagreement with the findings of Carnochan 
(1986), FM88 and JG93, the latter of whom found a correlation coefficient r— 0.33 for these 
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parameters. These groups suggested that wider bumps tended to be found along sighthnes 
with steep FUV rises, though this correlation was weak. In the present sample, there was 
no correlation, regardless of whether the parameter C4 was normalized to E(B-V) or Ay. 
This reinforces the idea that factors influencing the bump width are distinct from those 
influencing the carrier of the FUV rise (FM88). 

Correlations between C2, C3, C4, and Rv~^, noted by CCM, are confirmed here, as is the 
absence of a correlation between xq and Ky^^. Also, while 7 shows real variation, it does 
not appear to be linked with Ry~^. These can be seen in Fig. 4. 

3.2. The 2175 A Bump 

As the 2175 A bump is the only spectral feature yet known in the UV, parameters which 
describe it were carefully considered. In Fig. 5, a relationship between Cs/Ry and 7, very 
similar to the well-known relationship between C3 and 7, can be seen. The values of Ca/Ry 
are more confined at lower 7, but as 7 increases, Ca/Ry widens its range, with a general trend 
to increase with 7. This was noted by both FM88 and JG93 in their studies of C3 and 7; 
JG93 attributed it to the fitting procedure, but FM88 suggested that these two parameters 
are truly related in some way. 

The width of the bump shows real variation and environmental dependence. Values for 
7 ranged from 0.63±0.03 /xm"^ (HD 24263) to 1.47±0.05 /^m'^ (HD 29647). This is a wider 
range than that reported by FM86 (7 =0.77±0.09 - 1.25±0.07 /xm"^), perhaps reflecting 
the larger volume of the Galaxy and wider range of environments covered by this study. 
The average value in our sample is 7 =0.92±0.12 pim-^. HD 29647 and HD 62542 have 
the broadest bumps, with 7=1.467±0.049 //m"^ and 7 = 1.304 ± 0.04 /xm~^, respectively. 
The average 7 for dense and diffuse sightlines as defined above do not differ significantly, as 
7dense ~ ^-^^ ±0.04 ^in~^ and 7^j'Jj„se = 0.87 ±0.03 fim~^. This likely reficcts the density 
parameter's inability to sort out high versus low densities over long distances, as discussed 
previously. A z-tcst (i.e., Naiman, Roscnfeld, & Zirkel 1983) was done on the diffuse and 
dense 7 subsets. A z-test is essentiaUy the same as a t-test, as it is used to determine if 
two different datasets have significantly different mean values, but is used for large samples. 
As with the t-test, the resulting significance of such a test ranges from (the samples 
have significantly different means) to 1 (the samples have essentially the same means). For 
example, a significance of a few hundredths indicates that two samples have significantly 
different means. For the dense and diffuse sightlines, the resulting significance was ~ 10~^, 
which implies that there is a very significant difference between the mean 7 values of the 
dense and diffuse subsets. Fig. 6 shows a plot of all of the sightlines and their environmental 
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dependence, and a plot of sightlines which are more than 3a from the mean value of 7. It 

can be seen that there is a clear trend that the narrowest bumps tend to be found along 
sightlines with Ay/d <0.9 mag/kpc, the "diffuse" subset. These results strengthen those 
of FM86 and Cardelli & Clayton (1991), who found that lines of sight that passed through 
bright nebulosities ( "diffuse" sightlines) had narrower bumps than those that passed through 
dark clouds. 

Fig. 7 compares C4/R(V) and C3/R(V), with respect to 7. It can be seen that as the 
bump height increases, the range of values of the FUV rise increases, and the sightlines with 
the widest bumps tend to have higher values of C3/R(V), all of which has been reported 
previously by FM88. This is likely a reflection of the dependence of C3/R(V) and 7 on R^^. 
No similar trends can be seen between C4/R(V) and 7, as was discussed earlier. 

Fig. 8 plots 7 vs. xq, the bump central wavelength, for the sample. This figure shows 
that while 7 varies significantly from one sightline to another, xq does not. The extremely 
narrow range of xq agrees with FM86, who found that the mean xq = 4.60 //m~^. They 
find extreme values of xq are within ±0.04 fj,m~^ of the mean. Neither HD 62542 nor HD 
29647 were included in the FM86 sample. JG93 found a similar value of <xo >= 4.58 ±0.01 
/xm~^ for the Aiello et al. (1988) sample. For the much larger sample studied here, <xo > 
= 4.59±0.01 fj,m-\ with values ranging from 4.50±0.04 fxm.-^ (HD 145792) to 4.70±0.03 
//m-i (HD 29647). 

With the possible exception of HD 29647, it appears that xq does not vary from sightline 
to sighthne. CardeUi & Savage (1988) found that the 2175 A bumps toward both HD 62542 

and HD 29647 were significantly shifted to shorter wavelengths. However, we find xq = 4.54± 
0.03 fim.-^ for HD 62542. Cardelh & Savage (1988) found xo= 4.74±0.03 /im^^ and in the 
extinction curve shown in their work, the bump is visibly shifted blueward. In an effort to 
reproduce their results, the same lUE spectra and UV comparison star were used to construct 
an extinction curve. This curve was then fit over three ranges; the original FM90 range (3.3 
- 8.7 /xm^^), from 3.7 - 8.0 /xm~^, and from 3.7 - 8.7 /xm~^ excluding the region around Lya. 
None of these ranges produced a shifted bump. While the spectra used by Cardelli & Savage 
(1988) are the same as used here, their spectra were reduced and calibrated with a different 
software package than those in the final archive. 

In our sample of 417 sightlines, only one was shifted as much as 3 cr beyond the mean. 
We find xo = 4.70±0.03 //m"^ for HD 29647, the same value found by CardeUi & Savage 
(1988). The other two fitting ranges that were considered when testing for HD 62542's bump 
shift were applied to HD 29647. The values of xq that were found were within 1 a of the 
mean. Thus, it is quite possible that the shift in the HD 29647 bump also is fit dependent. 
Another sightline, lying close to that of HD 29647 has been studied (Clayton et al. 2003). 
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Both sightlines pass through the Taurus Dark Cloud. The second sighthne, toward HD 
283809, shows no evidence of a shifted bump. 

If the carrier is small graphite grains, then 7 oc 1/a, where a is the grain radius, while 
xq is not a function of size (Bohren & Huffman 1983; Hecht 1986). Hecht concluded that 
some fraction of small (a < 50A) carbonaceous grains are bare and these are responsible 
for the bump. In this scenario, the bump width is dependent on the temperature and 
size distribution of these grains. The remaining small carbon grains are hydrogenated, 
which would suppress the bump (Hecht 1986). One drawback to this model was that it 
required an overabundance of small grains in dense regions relative to diffuse regions (FM86; 
Sorrell 1990), since a large value of 7 implies small grain size in this scenario. This conflicts 
with the observation that dense regions tend to have higher values of Ry, indicating the 
prevalence of large grains. However, since xq is not dependent on Ry (or anything else), 
the grains which produce the bump may form a separate population from those which are 
responsible for variations in Ry (CCM). On the other hand, it can be seen in Fig. 4 that 
there is a weak correlation between 7 and Ry^^. This agrees with the finding that 7 is 
environment dependent. In 1990, Sorrell expanded upon Hecht's (1986) work, suggesting 
that the accretion of hydrogen on larger graphite grains with 60 A< a < 80 A was responsible 
for the variations in 7. He suggested that hydrogenation of grains in dense regions broadens 
the bump, without affecting xq. 

Mathis (1994) considered diamond, amorphous carbon (AMC), water ice, and PAH 
mantles on graphite cores. He found that the diamond coating tended to shift xq to longer 
much wavelengths. A thin coating of amorphous carbon yields xq =4.57 /^m"^ but produces 
a low 7= 0.84 ij,m~^. Increasing the thickness of mantle does not increase 7 to values that 
are observed; a thicker mantle only broadens the bump to 7=0.93 iim~^. Hydrogenated 
amorphous carbon mantles are even less suitable, as they shift xq to longer wavelengths 
while having even less of an impact on 7 than AMC did. Water ice mantles may broaden 
7, but only one sightline in our sample (HD 29647) has the telltale 3.07 fim ice feature 
(Goebel 1983). Warren (1984) showed that ice mantles do not affect the bump because there 
is almost no absorption at wavelengths near 2175 A. Neutral PAHs are also possible bump 
grain candidates, as they can produce both xq and 7 consistent with observations, with xq ~ 
4.61 /xm~^ for 7= 1.0 /xm~^ (Mathis 1994). However, these values were found by assuming 
that the optical constants of isolated PAHs are similar to those comprising the mantle and 
by ignoring impurities which might alter the PAH optical properties. Another drawback is 
that neutral PAHs should have absorption features in the UV which are not seen (Clayton 
et al. 2003), though Joblin et al. (1992) showed that this could be masked by averaging over 
a distribution of PAHs of different sizes. Also, it is unlikely that many PAHs are neutral in 
the diffuse ISM (e.g., LePage et al. 2003). 
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3.3. Deviations from CCM 

Adherence to the Ry-dependent extinction law of CCM in our sample was also examined. 
In order to do this, we found the CCM curve which best fit each extinction curve. The 
extinction curves were fit with a CCM curve using a standard IDL curve fitting routine that 
minimized the . The values of Ry found this way were compared to those calculated using 
IR photometry as is shown in Fig. 9. While there is some scatter, there is generally good 
agreement between the two values, with the fitted Ry being within Scr of the calculated 
value for 93% of the sample. This indicates that the vast majority of sightlines agree with 
CCM within the measurement uncertainties. The near absence of non-CCM sightlines in 
the Galaxy could be attributed to the lUE dataset somehow favoring lines of sight through 
CCM-type dust environments. lUE was a small telescope with limited dynamic range. As 
a result, there is little information on grain parameters deep inside dark clouds. The dust 
columns sampled in this study all have Ay < 5 mag. 

Next we measured deviations of each measured extinction curve from its best fit CCM 
curve. We followed the method of Mathis & Cardelli (1992; hereafter MC92), in which the 
authors specified the deviation at various wavelengths to be: 

b{xi) = [Ax/Av]i - [Ax/Av]ccM 

Plots are shown of these deviations in Fig. 10 at 4.65, 4.90, 5.07, 5.24, and 7.82 /i,m~^. The 
calculated and best-fit CCM extinction curves agreed within 2a of each other for all but four 
of the sightlines, HD 210121, HD 204827, HD 62542, and HD 29647, which are deviant at 
several wavelengths. Their deviations are shown in Table 7. 

HD 204827 is an intriguing sightline. It has been shown that when a foreground red- 
dening component has been removed, its extinction curve resembles that of the SMC, even 
though the environment of the dust in the HD 204827 cloud is quite different from that 
seen in the SMC sightlines (Valencic et al. 2003). A similar curve has been found toward 
HD 283809 (Whittet et al. 2004). The dust local to HD 204827 has high density. The HD 
204827 dust cloud resembles a molecular cloud more than the diffuse ISM. The HD 204827 
cloud is very rich in carbon molecules, showing large column densities of C2, C3, CH, and 
CN (Oka et al. 2003; Thorburn et al. 2003). The sightline to HD 204827 samples dust swept 
up by a supernova or hot star winds (Patel 1998). In these respects, the HD 204827 dust 
is quite similar to the sighthne toward HD 62542 (CardeUi & Savage 1988). Its dust is also 
rich in carbon molecules. The sighthne to HD 62542 also lies on the edge of material swept 
up by a stellar wind bubble. The other two non-CCM, weak bump, steep far-UV sightlines 
in the Galaxy, HD 29647, and HD 210121, are also associated with dense clouds (Cardelli 
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fe Savage 1988; Larson, Whittet, & Hough 1996; Cardelli & Wallerstein 1989). The dust in 
the molecular cloud associated with HD 210121 is likely to have been processed as it was 
propelled into the halo during a Galactic fountain or other event. The final sightline, toward 
HD 29647, is sampling dust in a quiescent dense cloud. The steep far-UV extinction in these 
clouds helps shield the molecules in these clouds from dissociating UV radiation leading to 
larger column densities (Mathis 1990). 

Fig. 11 shows the deviations, 5(7.82) and 5(4.65), plotted against the column density 
of CN along the sightline divided by Ay for 46 sightlines from our sample, including the 
four non-CCM sightlines discussed above (Fcderman 1994; Oka ct al. 2003). CN can be 
used as an indicator of relatively dense regions in diffuse, molecule-rich clouds (Joseph et al. 
1986; Gredel et al. 2002). The average 5(7.82) and N(CN)/Ay (lO^^cm-^ mag"^) for the 
42 CCM-hke sighthnes were 0.03±0.03 and 0.33±0.01, respectively. For the four non-CCM 
sightlines, these values were 1.06±0.15 and 2.75±0.15. These four sighthnes have both high 
N(CN)/Ay and sigificant deviation from CCM at a:=7.82 compared to the CCM-like sample. 
Similarly, Burgh et al. (2000) showed that sightlines with steep FUV rises tend to have high 
N(CO)/E(B-V), though they did not notice any correlating weakness in the bump strength. 

Virtually all of the non-CCM sightlines known today are in the (diffuse) Magellanic 
Clouds, not in dense Galactic clouds (e.g., Gordon ct al. 2003). Clayton et al. (2000), found 
that the average of seven low density, low reddening sightlines showed an extinction curve 
very similar to the LMC. These 7 stars were found to be behind gas that showed signs of 
being subjected to shocks The deviations from CCM for the average of these 7 sightlines as 
well as for the SMC and LMC average curves were also found and are included in Figs. 10 
and 9. Similar deviations from CCM may arise in various dust environments. 



3.4. The Updated CCM Law 

The new extinction curves and values of Ky allow us to revise the CCM law. This was 
done by following the method described in CCM; that is, by plotting Ax/ Ay against Ry^ 
for all wavelengths, then performing a linear least-squares fit of the resulting correlation. 
Thus, for all wavelengths, the extinction Ax/ Ay = a{x) + b{x)/Ry. Then, a{x) and b{x) 
were plotted against x and fit with a polynomial expression by minimizing x^. It was found 
that for 3.3 > X > 8.0 //m~^. 



a{x) = 1.808 - 0.215X - 



0.134 



(x- 4.558)2 + 0.566 



+ Faix) 
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b{x) 



= -2.350 + 1.403X + 



1.103 



{x - 4.587)2 + 0.263 



where, for x < 5.9/im 



1 



Fa{x) = Fb{x) = 0.0 



and for 5.9 < x < 8.0/im 



1 



Fa{x) -- 



-0.0077(x - 5.9)^ - 0.0003(x - 5.9)=^ 
0.2060(x - 5.9)^ - 0.0550(x - 5.9)^ 



A comparison of extinction laws are plotted in Fig. 12. The most noticeable difference 
between the law found here and the original CCM law, the increase in the zero-point of the 
UV extinction, is a consequence of the slightly higher value of the first term in a{x). In 1988, 
CCM also found a higher value for this term (1.802) than in 1989 (1.752), when they found 
that the lower value, which was within the dispersion in the data, was needed for the curve 
to smoothly join the optical extinction data. Fig. 12 also shows the resulting curve when the 
average FM parameter values are used to construct extinction curves at different values of 
Ry, in a fashion similar to Fitzpatrick (1999). The two methods of curve construction yield 
curves that are within 5% of each other for 2.5 < Ry < 5.0 over the entire wavelength range. 
For all values of Ry, the difference is primarily in the bump height, or the FM parameter 
C3. This is most noticeable for 5.0 < Ry < 6.0; here, the curves agree to within 7%, with 
the curve found through the average FM parameters having a weaker bump than predicted 
by the either the revised or original CCM law. This can also be seen to a certain extent in 
Fig. 4. A weakened bump at high Ry may be expected from Whittet et al.'s (2004) work, 
which suggests that the bump carrier may be destroyed in dark clouds. Both the revised 
CCM curves and the average FM curves are within about 5% of the original CCM law for 
2.5 < Ry < 3.5. However, for increasing values of Ry, the difference between the curves 
increases as well, reaching 20% at Ry ~ 5.0. That the original CCM law agreed so well with 
the extinction curves examined here is an indicator of the sizes of the uncertainties on the 
data. Changes in the value of the extinction zero-point is an important factor. If the first 
term in a{x) is decreased by 0.20, the new CCM law and the original curve agree to within 
15% for 2.5 < Ry < 6.0 for x < 7.0 nm'^, and within about 20% for x > 7.4 iJ,m-\ The 
dotted curve, constructed from Fitzpatrick's (1999) suggested values of FM parameters with 
only C2 (and through it, Ci) dependent on Ry, illustrate further the dependences of some 
FM parameters on Ry. 
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4. Conclusions 

We have constructed a homogeneous database of UV extinction curves, all of which 
have been fit by the FM relation and for which values of Ry have been calculated. All 
extinction curves and FM parameters have been normalized to Ay, rather than E(B-V), 
so that relationships between parameters may be more easily seen. This is the largest and 
most comprehensive database of extinction curves yet constructed, sampling a wide variety 
of environments. It contains much longer sightlines than those included in previous studies 
(30 sightlines have d > 5 kpc), thus greatly increasing the volume of the Galaxy sampled 
compared to previous works. Over 150 sightlines have d > 2 kpc, so the regions of the Galaxy 
beyond the solar neighborhood are well represented. 

The main results of this study are: 

(1.) The GGM extinction relation accurately describes the diffuse Galactic ISM in virtually 
all cases. Out of 417 sighthnes, only 4 deviated consistently from their GGM extinction 
curves by more than Sa. This indicates that the physical processes which give rise to grain 
populations that have CCM-like extinction dominate the ISM, and thus, the quantity Ry 
can accurately describe the UV extinction for most sightlines. This implies that the grain 
populations responsible for different components of the extinction curve are being processed 
efficiently and systematically along most sightlines. 

The new curves and Ry values allow for an updated GGM law to be made. The original 
and updated versions are within 5% of each other for 2.5 < Ry < 3.5 over the full wavelength 
range covered, 3.3 < x < 8.0 /xm"^. At higher values of Ry, the differences in the extinction 
zero-point reduce the agreement in the curves to about 20%. 

(2.) The bump width has a strong environmental dependence, with narrow bumps favoring 
diffuse sighthnes, and broad bumps favoring dense sighthnes. Very broad bumps are rare, 
as only 21 out of 417 sighthnes had 7 > 1.1. 

(3.) The central wavelength of the bump is invariant and may be regarded as a constant, with 
xo=4.59±0.01. Unhke the other parts of the UV extinction curve, this parameter does not 
respond to different environments. The invariance of xq, along with the observed variations 
in 7, put strong constraints on possible bump grain and mantle materials. 

(4.) While there is evidence for shock processing in three of the four non-GGM sightlines, the 
common denominator is that all four sightlines have dense, molecule-rich clouds. They also 
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had weak bumps and strong FUV extinction for their measured Ry CCM curves, especially 
those which pass through cold, quiescent regions. The weakened bump may reflect processing 
which modifles or destroys the bump carrier in dark clouds (Whittet et al. 2004). This is 
seen in the bumpless dust toward HD 204827 and HD 283809. Other dense sightlines with 
strong FUV extinction, HD 204827, HD 62542, and HD 210121, pass through dense clouds 
that may have been exposed to shocks or strong UV radiation that disrupt large grains, 
possibly resulting in a size distribution that is skewed toward small grains. 

Together, the sightlines in the Galaxy and the Magellanic Clouds suggest that similar 
extinction properties may arise from very different environments. 

Financial support was provided by the Louisiana Board of Regents, BoRSF, under agree- 
ment NASA/LEQSF(1996-2001)-LaSPACE-01 or NASA/LEQSF (2001-2005)-LaSPACE and 
NASA/LaSPACE under grant NGT5-40115. 
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Table 1: Comparison of Previous UV Extinction Studies to the Present Work 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Authorship 


No. of Sighthnes 


Instrument 


FM 


Spectral 


Ry 










Classification 




Savage et al. 1985 


1415 


ANS 


no 


MK 


no 


Aiello et al. 1988 


115 


lUE 


no 


UV 


no 


FM90 


78 


lUE 


yes 


UV 


no 


Papaj et al. 1991 


166 


TD-1 


no 


UV 


no 


JG93 


115 


lUE 


yes 


UV 


no 


Barbaro et al. 2001 


252 


ANS 


yes 


MK 


yes 


This Work 


417 


lUE 


yes 


UV 


yes 



Table 2: Photometry and Sources Used in Present Work 



Name 


Spectral 


uv 


U 


B 


V 


J 


H 


K 


UBV 


U Non-2MASS 




Type 


Type 














Source 


Source^ Source 


HD 14052 


Bllb 


Bllab 


7.90 


8.48 


8.18 


7.51 


7.50 


7.42 


1 










±0.01 


±0.01 


±0.01 


±0.03 


±0.03 


±0.02 






HD 14357 


B2II 


B2Ib 


8.31 


8.83 


8.52 


7.86 


7.83 


7.80 


1 










±0.01 


±0.01 


±0.01 


±0.03 


±0.02 


±0.02 






HD 41690 


BIV 


B1.5in 


7.31 


7.93 


7.71 


7.20 


7.23 


7.20 


1 










±0.01 


±0.01 


±0.01 


±0.02 


±0.02 


±0.03 






HD 199216 


Bin 


B2III 


7.05 


7.51 


7.03 


6.05 


6.03 


5.97 


1 










±0.01 


±0.01 


±0.01 


±0.02 


±0.03 


±0.02 






HD 203938 


B0.5IV 


B1.5IH 


7.13 


7.54 


7.08 


5.97 


5.88 


5.81 


1,2 










±0.01 


±0.01 


±0.01 


±0.02 


±0.02 


±0.02 






HD 235783 


Bllb 


Bllab 


8.14 


8.85 


8.68 


8.28 


8.27 


8.26 


1 










±0.01 


±0.01 


±0.01 


±0.03 


±0.05 


±0.02 







' If not the same as UBV Source. 

The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. 



Table 3: Photometry Sources 



Number Citation 

1 Hiltner, W. A. 1956, ApJS, 2, 389 

2 Nicolet, B. 1978, A&AS, 34, 1 

3 Schild, R., Garrison, R., & Hiltner, W. 1983, ApJS, 51, 321 

4 Haupt, H., & Schroll, A. 1974, A&AS, 15, 311 

5 Guetter, H. 1974, PASP, 86, 795 



The complete version of this table is in the electronic edition of the Journal. The printed 
edition contains only a sample. 



Table 4: General Extinction Characteristics 



Name 


E(B-V) 


Ry 


Ay 


Distance 


Max. Distance 


Min. Distance 










(pc) 


(pc) 


(pc) 


HD 14052 


0.49 


2.99 


1.47 


4195.08 


3322.81 


6065.88 




±0.04 


±0.19 


±0.14 








HD 14357 


0.48 


2.90 


1.39 


3817.16 


2838.24 


5229.96 




±0.04 


±0.21 


±0.15 








HD 41690 


0.47 


2.91 


1.37 


820.92 


602.25 


1111.31 




±0.04 


±0.22 


±0.16 








HD 199216 


0.72 


2.62 


1.88 


1567.69 


1155.49 


2380.10 




±0.04 


±0.36 


±0.27 








HD 203938 


0.71 


3.13 


2.22 


496.81 


362.62 


672.56 




±0.04 


±0.14 


±0.16 









The complete version of this table is in the electronic edition of the Journal. The printed 
edition contains only a sample. 
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Table 5: FM Parameters 



Name 


ci/Ry + 1.0 


C2/Ry 




C4/Ry 


Xo 


7 


HD 14052 


1.1856 


0.2271 


0.8595 


0.2124 


4.5740 


0.8440 




±0.3918 


±0.0258 


±0.1237 


±0.0411 


±0.0140 


±0.0280 


HD 14357 


1.1134 


0.2745 


1.0155 


0.2421 


4.5830 


0.8970 




±0.1461 


±0.0282 


±0.1364 


±0.0296 


±0.0050 


±0.0280 


HD 41690 


1.3412 


0.2326 


0.7100 


0.2790 


4.6050 


0.7800 




±0.2734 


±0.0347 


±0.1193 


±0.0494 


±0.0120 


±0.0270 


HD 199216 


0.7992 


0.2908 


1.3424 


0.1408 


4.5900 


0.9360 




±0.2038 


±0.0436 


±0.2244 


±0.0339 


±0.0090 


±0.0340 


HD 203938 


1.1581 


0.1933 


1.1735 


0.1371 


4.5580 


1.0120 




±0.3085 


±0.0190 


±0.1534 


±0.0217 


±0.0110 


±0.0340 



The complete version of this table is in the electronic edition of the Journal. The printed 
edition contains only a sample. 



Table 6: Correlations between FM Parameters 



r 


Absolute Deviation 


Linear Fit 


0.58 


0.19 


C3 = 2.12 7 - 0.95 


-0.79 


0.13 


ci = -2.54 C2 + 1.67 


0.49 


0.54 


C3 = 1.75 C2 ± 0.54 



Table 7: Deviations of Non-CCM Sighthnes 



HD 


5(4.65) 


5(4.90) 


5(5.07) 


5(5.24) 


5(7.82) 


29647 


-0.58 


-0.35 


-0.18 


-0.08 


0.31 


62542 


-1.49 


-1.10 


-0.75 


-0.50 


1.12 


204827 


-0.65 


-0.44 


-0.29 


-0.18 


0.43 


210121 


-2.63 


-2.35 


-1.76 


-1.21 


2.33 
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Fig. 1. — Top panel: All siglitlines with published FM parameters. Diamonds: Aiello et al. 
(1988); Triangles: FM90; Squares: OB associations with at least one member having FM 
parameters. Bottom panel: All sightlines in the database. Diamonds: field stars. Squares: 
OB associations with at least one member having FM parameters. The Galaxy overlay is 
from Vallee (2002). 




Fig. 2. — FM parameters of database sightlines. All parameters, except xq and 7, have been 
divided by Ry and thus normalized to Ay- Dense sightlines (those with Ay /d> 0.9mag/kpc) 
are indicated by triangles, while diffuse sightlines {Ay/d < O.Qmag/kpc) are indicated with 
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Fig. 3. — Top panel: The FM parameters C2 and C3 from all database sightlines, divided by 
Ry and thus normalized to Ay. Bottom panel: Values of C2 and C3 from FM88 (triangles) 
and JG93 (squares), divided by Ry and thus normalized to Ay. Representative error bars 
are indicated in both panels. 
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Fig. 4. — FM parameters plotted against Rv~^- All parameters, except xq and 7, have been 
divided by Ry and thus normalized to Ay- Symbols are the same as defined in Fig. 2. The 
solid lines represent the best fits. The dashed lines indicate the expected values from the 
CCM relation; the dash-dot lines are from Fitzpatrick's (1999) reformulation of the CCM 
law. Representative error bars are indicated. 
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Fig. 5. — Bump height, C3, vs. bump width, 7. Symbols are the same as defined in Fig. 2. 
Representative error bars are indicated. 




Fig. 6. — Bump width, 7, vs. density. The soHd hnc is the cutofi^ between dense and difi^use 
fines of sight. Left panel: all sightlincs in database. Right panel: Sightlines beyond Scr of tfie 
mean. In order to make the correlation easier to see, those lines of sight that agreed with 
the mean were removed; the remainder are shown in the right hand panel. In both, the line 
at Ay/d =0.9 indicates the cutoff between dense and diffuse lines of sight. 
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0.4 - 




-0.2 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 
c/R(V) 



Fig. 7. — FUV curvature, C4/Ry, vs. bump height, Ca/Ry. Symbols are the same as defined 
in Fig. 2. Representative error bars are indicated. The bump height (ca/Ry) and FUV 
curvature (c4/Ry) are plotted with respect to bump width. Values of 7 were split into 
three categories, each roughly pertaining to environment, with diamonds indicating 7 < 0.9, 
triangles for 0.9 < 7 < 1.1, and squares for 7 > 1.1. 



4.90 




0.6 0.8 1.0 1.2 1.4 1.6 

? 



Fig. 8. — Bump central wavelength, xq, vs. bump width, 7. The solid line is the average. 
Symbols are the same as defined in Fig. 2. Representative error bars are indicated. 
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Fig. 9. — The best-fit Ry values, found through minimization, compared to Ry vahies 
found through IR photometry. A Une of slope unity has been overlayed. Four lines of sight 
with unusual extinction are indicated (filled circles). Lines of sight where the best-fit Ry 
did not agree with the measured value are indicated with error bars. The average LMC 
(filled triangle), LMC2 Supershell (filled diamond), and SMC (filled star) are plotted. Prom 
Clayton et al. (2000), the average low-density sightline for a Galactic region with extinction 
similar to that of the Magellanic Clouds (the "SD region" ) is also shown (filled square) . 



-30- 



-2 
6(4,65) 



-2 
6(4.65) 



-2 
6(4,65) 



-2 
(5(4,65) 



-2,0 -l,b -1,0 -0,b 0.0 0.5 1.0 
6(5.24) 



Fig. 10. — The deviations from Ry law at various wavelengths. Representative error bars are 
indicated. Four Galactic lines of sight with unusual extinction are named. The SD region 
average sightline (Clayton et al. 2000; filled square), average LMC (filled triangle), LMC2 
Supershell (filled diamond), and SMC (filled star) are also plotted. 
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Fig. 11. — Deviation at 7.82 /xm~^ and 4.65 fira^^ versus abundance of CN. Values of N(CN) 
are from Federman (1994) and Oka et al. (2003). Unusual Galactic sightlines are indicated 
with solid cirles. 
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Fig. 12. — A comparison of extinction laws. Solid line: the extinction law derived in this 
work; dashed line: the original CCM law; dash-dot line: the curve constructed with the 
FM parameter averages found in this work; dotted line: the curve constructed with the FM 
parameters suggested by Fitzpatrick 1999. 



